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Abstract. The Dulac criterion is a classical metliod to rule out the existence 
of periodic solutions in planar differential equations. In this paper the applica- 
bility and therefore reversibility of this criterion is under consideration. 



o 

1. Introduction AND Motivation 

^— • 

O Let X : V ^ E? he a a smooth vector- valued function with components P 

w and Q, X = {P,Q), defined on some planar domain D C jsjow consider the 

Q\ system of two ordinary differential equations 

(1-1) ^ = p{x,y), f = Q(^,y)- 

Q Such systems appear frequently in applications, for example in electrical engi- 

^ neering, physics, biology and many others, but they have also become a field of 

^ mathematical interest on their own. For simplicity, one can think of P and Q either 

^ just as polynomials of two real variables x and y, or as smooth or even analytical 

I— I functions given by some power series that converges in V, i.e., X is of class C 

^ with r G {1, 2, . . . oo,uj}. Geometrically speaking, the solutions of system ( |1.1[ ) 

^ are smooth curves in the plane that are tangential to the function X at each point 

Un z = (x, y) G V. If the components of X are polynomials in x and y of degree one 

^ with a,b,c,d € M. being real constants, the analysis of system ( |1.1[ ) is easy and 

the solution curves of system ( |1.1| are given explicitly in terms of the matrix ex- 
y—i ponential function e"^ := YlT=o M^'^- "^^^ situation changes drastically if X is 

^ at least quadratic, i.e., if P and/or Q are polynomials of degree two or higher, 

P,Q& M.2[x,y]. Many things are known about quadratic differential equations, 
see |[T][TT}[T2l, but they are still a broad area of research. Quite a few people dedi- 
cated their whole life on the investigation of quadratic planar differential equations. 
The reason for that is related to the 16-th Hilbert problem, which is still unsolved 
even in the quadratic case. 

1.1. Sixteenth Hilbert problem. Let Rd[x,y] denote the space given by poly- 
nomials in X and y of degree at most d and consider a planar system ( |1.1| ) with 
polynomial right-hand-side X G ]R^[x, y]. The 16-the Hilbert problem consists of 
two parts, the first is a purely algebraic problem and deals with questions related 
to the topology of algebraic curves and surfaces. In this paper we are mainly con- 
cerned about the second part of the 16-th Hilbert problem which deals with the 
curves that arise as solutions of the planar differential equation ( |1.1| ) with poly- 
nomial vector field X. More precisely, in the second part of his problem Hilbert 
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asks for the maximal number and relative position of the isolated closed orbits, 
called limit cycles the differential system ( |1.1| ) can have at most. It is even diffi- 
cult to understand why a fixed planar differential equations with polynomial vector 
field can only have a finite number of limit cycles. A proof of this fact is due to 
Ilyashenko in 1991 \ who actually corrected a very complicated and long proof 
due to Dulac a student of Poincare. Ecalle et al. obtained parallel and indepen- 
dently a proof for this theorem [4]. Note that this does not imply the existence of an 
upper bound for the maximal number of limit cycles H{d) a planar polynomial sys- 
tem of degree d can have. One calls H{d) the d-th Hilbert number. Linear vector 
fields have no limit cycles; hence H{\) = 0. Quadratic systems can actually have 
four limit cycles and some people believe that this is the maximal number of limit 
cycles a quadratic system can have, but it is still unknown whether or not H{2) 
is a finite number. Usually, the first part of the 16-th Hilbert problem is studied 
by researchers in real algebraic geometry, while the second part is considered by 
mathematicians working in dynamical systems or differential equations. Hilbert 
also pointed out that there exist possibly connections between these two parts. See 
|[8| and Q for the original paper in Russian and more recent | [T3| for a survey about 
the second part of the 16-th Hilbert problem. 

1 .2. Dulac criterion. In 2008 the author was confronted with the analysis of some 
quadratic polynomial differential equation. In literature there existed a quite techni- 
cal proof showing the uniqueness of a limit cycle for this system. Then, fortunately, 
he was able to obtain the same result by applying the well-known 

Theorem 1 {Dulac criterion). Let Q. (1 ^ be a simply connected region, X := 
{P,Q) S C^(ri,M^) a smooth vector field and B £ C^(r2,M) a smooth real-valued 
function - such that the partial differential inequality 

(1.2) d„(BX):='^ + «(|<5i>0 

OX oy 

is satisfied in ^1. Then the planar ordinary differential equation ( |1.1[ ) with X as 
right-hand-side does not posses any periodic solution that is fully contained in Q. 

Remark 1. (1) The proof of theorem^is indirect: One assumes the existence 
of a closed orbit of ( |1.1| ) in and applies the divergence theorem of GauB. 
This gives immediately a contradiction. 

(2) In 2009 the autorwas able to weaken the assumptions made on the function 
B, see |7|. Essentially, the statement of theorem^still holds if the function 
B is only weakly differentiable of first degree and equation ( |1.2| ) holds only 
almost everywhere in 0,, compare definition^ 

(3) If the domain C referred to in theorern^is not simply connected, but 
p-connected for some p € N, p > 2, then there can be at most p—1 closed 
orbits fully contained in Q. 

(4) The Dulac criterion is a generalization of the Bendixson criterion. In fact, 
the Bendixson criterion follows from theorem^if we choose B = 1, the 
function that is constant Ifor all z £ Q. 



^During history the term hmit cycle was used for a stable isolated closed orbit or just for closed 
orbit. Nowadays and in this paper, a limit cycle is a closed orbit that is isolated in the set of all 
periodic orbits of some differential equation. 
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Definition 1 (Dulac function). The multiplier B € Ty^'^(J7,M) is called Dulac 
function of the planar dynamical system \\.\) in Q D if and only if there is a 
real-valued continuous function g > having positive sign in Q except on a set of 
Lebesgue measure zero such that equation 

(1.3) div(SX) = 5-divX + (V-B,X) =5 

holds in Q.. Here VF^'P(f],M) denotes the Sobolev space of LP functions that are 
weak differentiable or first degree. Vice versa, we call any function of class W^'^ 
satisfying ( |1.3| ) a g-DxAac function of X in il. 

In general, it is very difficult to find a Dulac for some given vector field X, but 
if one incidentally guesses such a function, it is quite easy to verify that it obeys 
condition ( |1.2| p| Just for curiosity, the author wondered whether the converse state- 
ment of theorem[T]is also true. 

Question 1. Given a smooth planar vector field X and assume that the planar 
differential equation ( |1.1[ ) does not have any periodic solution in some simply con- 
nected domain Q.. 

Does there exist a Dulac function in Q? 

The answere is No if the boundary of is formed by a periodic solution of the 
corresponding system, but in every other case that was considered the author was 
able to obtain an affirmative answere to question [T] from above. In the following 
we will to quote some of these positive results. 

2. Local results 

Gradient systems never possess periodic solutions, see |6|, hence we expect that 
gradient fields have a Dulac function. Indeed, this is the case. 

Theorem 2 (Gradient fields). Let X = W be a globally defined gradient field 
with potential V G C2(i^2^ ]^^)_ j^^^^ ^ ^y^^ ^ (^_y^ S3 = 1^ 

are local Dulac functions for X in some Qi C M^, i = 1, 2, 3, such that U^^^ilj = 

Proof. Define Bi = exp {V), B2 = exp {—V) and B^ = V and set ilj := {z E 
M^l div{BiX){z) > 0}. Verify that uf^^ili = gives indeed the whole plane. 

□ 

A standard result from calculus says that vector fields can be straightened locally 
unless there is an equilibrium. This fact can be used to obtain a local existence 
result in domains not containing zeros of the vector field. Such domains are called 
canonical regions. 

Theorem 3 (Parallel flow). Let X : V C M? ^ M? be any smooth vector field 
and Q, a canonical region of X, i. e. a region where the fiow of system ( |1.1[ ) is 
equivalent to a parallel flow in the sense of Neumann |T5] . Then, under some 
integrabilty assumptions, X has a Dulac function in Vl. 



One can compare it with the decomposition of a natural number into its prime factors: If some 
natural number n together with some product p :— HiLi Pi^ of powers of prime numbers is given, 
it is quite easy to decide whether p is the prime factorization ofn,p = n, but, from our current state 
of knowledge, it is very difficult to decompose a big natural number into its prime factors. 
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Proof. By assumption there are no equilibria in the domain, hence the vector field 
can be straightened locally. Observe that, by this process of straightening, the 
planar system ( |1.1[ ) decouples and the partial differential equation ( L3 1 from the 



definition of a Dulac function reduces to the linear ordinary differential equation 

BdiYvX + aBr = g 

which implies 

g-BdivX 

Br = 

a 

This is an ordinary differential equation with solution 

where a{r) = J div X dr and k G R is a constant of integration. □ 

For linear vector fields there is also a Dulac function in terms of a quadratic 
function if some spectral condition holds. 

Theorem 4 (Linear case). ||7| Let X{z) = Az be a linear vector field with A 
having at most one eigenvalue zero. The linear system ( |1.1[ ) does not have periodic 
orbits if and only if the spectrum of the matrix A consists only of eigenvalues with 
nonzero realpart or zero, i. e. (y{A) n iM C {0}. 

Proof. Let the vector field be X{z) = Az with 

c aj 

Instead of the general partial differential inequality ( |1.2| ) the author considered 

(2.1) div(SX) = := p2 _^q2 

and made a quadratic ansatz for the Dulac function 

B=\{z,Gz) 

with some matrix G G R^^^. Then equation ( |2.1[ ) reduces to 

A^G + GA + iTA-G = A^A. 
Letting S := A + tv A, one obtains 

(2.2) S^G + GS = A^A. 



Equation ( |2.2[ ) is the well-known Lyapunov equation, a special case of the more 
general Silvester equation. Now one can apply the well-established solvability 
theory for the Silvester equation and verify that the Lyapunov equation ( |2.2[ ) indeed 
has a unique solution under the spectral assumptions that were made. On the other 
hand, one obtains by direct calculation 

B = h2ox^ + 6o22/^ + hiixy + bi^x + boiy + boo 

with 

a'^ + 4a^d - a^{2bc - - 3d^) + 3acd{c - b) + c^{b^ -bc + d^) 



b20 



(o + d)(3a2 + lOod - 46c + 3^2) 



_ a2(b2 + 3rf2) ^ ad(362 - 36c + M"^) - b^c + b'^jc'^ + d^) _ 2bcd^ + 
°^ ~ [a + d)(3a2 + Wad - Abe + 3d'^) ' 
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_ 2a% + a^d{7b + 3c) - a{3b^c + b{c-3£) - led?) - cdjb'^ + 36c - 2£) 
" ~ (a + d)(3a2 + Wad - Abe + M^) ' 

and boi = bio = ^oo = 0. The case when tryl = a + d = has to be examined 
with care. The reader may verify that having spectrum o-{A) = {0, | tr ^4} is 

equivalent to det A = (^^)^. Hence the quadratic Dulac function from Q does 
the job because we assumed that at most one eigenvalue of A is zero. □ 

The Hartman Grobman theorem combined with the last two results gives some 
local existence statement holding in some neighborhood of hyperbolic fixed points. 

Proposition 1. Let X be a smooth vector field and z a hyperbolic zero of X, i. e., 
the realpart Re z ^ Q is different from zero. Then there is a neighborhood U of z 
such that X has a Dulac function on U. 



3. Qualitative theory 

Let us recall now some basic definitions and results that are frequently used in 
the qualitative theory of planar differential equations. For a more detailed intro- 
duction we refer to |[5| and |3| . We call system ( |1.1| ) integrable in the domain V 
if it has a first integral defined on this domain, i. e. a non-constant smooth scalar- 
valued function H of class which is constant on each solution (x(t), y{t)) of 
as long as it is defined. This means: if (x(t), y{t)) is any fixed solution of 
defined for t G [0, tm], its maximal interval of existence, and H G C^('D, M) 
a first integral of system ( |1.1| ), then there is a real number h such that 

(3.1) H{x{t),y{t)) = h for all t£lm 

is satisfied. Taking the derivative with respect to time t of equation p.l| ), we see 
that any first integral H ^ {V, M) of ( |1.1| ) satisfies the linear partial differential 
equation 

(3.2) {VH,X) =P-H^ + Q-Hy = 



Conversely, every non-constant solution of equation (3.2i is a first integral H 



P — M of ( |1.1| ). If Hq is a non-constant solution of p.2i, then every other solution 



is of the form F{Hq) where F is an arbitrary function having continuous partial 
derivatives (use the chain rule to verify this). First integrals are strongly related to 
the notion of integrating factors. 

Definition 2 (Integrating factor). An integrating factor /j, of the planar dynamical 
system is a smooth solution fj, G C^(il,R) of the linear partial differential 
equation 

(3.3) div(/i ■X) = fi-divX + {Vfi, X) = in Q. 



Note that the first equality in ( |3.3| ) is due to the Leibniz rule and div X, the di- 
vergence of the vector field X, denotes the trace of the Jacobian of X, div X = 

dP_ ,dQ 
dx ' dy ' 



If instead a first integral : D — )• M of \IA) is known, using equation p.2| ). 



one can equivalently define an integrating factor as the common value of the ratios 
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i. e., an integrating factor must satistfy both Hy = fiP and = —fJ-Q. The latter 
two relations can be read as 

(3.5) dH{x, y) = fi{x, y){P{x, y)dy - Q{x, y)dx), 

where dH{x, y) indicates the differential of H. Thus, multiplying the right hand 
side of by an integrating factor makes the equation an exact differential and 
an exact equation), which means that, whenever an integrating factor is available, 
the modified vector field jiX has vanishing divergence and the problem of solving 
equation ( |1.1| ) is reduced to one-dimensional integration 

H{x,y)= / fi{x,y){P{x,y)dy - Q{x,y)dx). 

Note that the latter line integral might not be well-defined if the domain V is not 
simply connected. For this reason, integrating factors are usually considered only 
in connected components of V. Secondly, we observe that the vector fields X and 
have the same phase portrait (with maybe reversed orientaion if // is negative) 
as long as /x does not vanish. That is why solving system ( |1.1[ ), i. e. constructing a 
first integral, and finding an integrating factor for it are considered to be equivalent 
problems. In applications the notion of inverse integrating factor is very common. 

Definition 3 (Inverse integrating factor). The function V G C^(il,M) is called an 
inverse integrating factor for the planar system ( |1.1[ ) in the domain E? if 

1 

is an integrating factor of system ( |1.1[ ) in \ {V = 0}. As usual {V = 0} is short 
notation for the preimage of zero under V, y^^(O) = {z G : V{z) = 0}. 

The method of integrating factors is, both historically and theoretically, a very 
important technique in the qualitative analysis of first order ordinary differential 
equations. The use of integrating factors goes back to Leonard Euler {1101 - 1783). 
Integration factors and first integrals refer to the problem of integrating the planar 
system ( |1.1| ), which means geometrically nothing but finding smooth curves that 
are tangential to the vector field at each point. Then, another interesting problem is 
deriving the qualitative behavior of these solution curves as, for instance, topology 
(whether they are closed or not), asymptotic (whether they blow up in finite time 
or remain in some compact set and approach a limit cycle) and stability. 

4. A UNIFYING POINT OF VIEW 

Definition 4 (Invariant curve). Let Q. be an open set of the plane. An invari- 
ant curve is the vanishing set or the primage of zero of some smooth function. More 
precisely, to a given function f £ C^($7, C) we associate the preimage of zero 

(4.1) /-i(0) = {/ = 0}:={zGl7|/(z) = 0} 

and call it an invariant curve of the vector field X = (P, Q) £ C^(2?, M^) if there 
is a smooth function k G C^($7, C), called cof actor of f, satisfying the relation 

(4.2) (V/, X) = k- f forallz^n 
or more explicitely 

f.,{z)-P{z) + fy{z)-Q{z) = k{z).f{z) 
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Here V/ = denotes the gradient of f, the function (•, •) is the canonical 

inner product of and the subscripts of f indicate partial derivates, fx = g^, 
fy = Note that on the invariant curve / = the gradient of f, V/, is 
orthogonal to the vector field X by the defining property of invariant curves ( |4.2| ). 
By convention, the function f defining the invariant curve /~^(0) C is called 
invariant function. 

For studying the multiplicity of invariant curves, the notion of exponential fac- 
tors, a special case of invariant functions, is useful. It allows the construction of 
first integrals for polynomial systems via the same method used by Darboux. 

Definition 5 (Exponential factor). Given two h,g^ M[x, y] coprime polynomials, 
the function exp{f^/g) is called an exponetial factor for system ( |1.1| ) if there is a 
polynomial k G M.[x, y] of degree at most d — 1, d := maxjdegP, degQ} being 
the degree of the polynomial system ( |1.1[ ), satisfying the relation 

(4.3) (V(e'/"),X) = /fc-e'/'^ 

Note that obviously {(e"/'") = 0} = 0, but g = defines an invariant curve for 
system ( |1.1| ). 

Definition 6 (Daboux function). Any function of the form 

r I 

i=i j=i 

where, for 1 < i < r and 1 < j < I, fi{z) = and gj{z) = are invariant curves 
for system \\A\ , hj is a polynomial of(C[x, y], Aj and fij are complex numbers and 
Tij is a natural number or zero, is called Darboux function. 

The following remark reminds the reader of some facts about invariant curves. 

Remark 2. (1) Invariant curves are very important in the qualitative study of 
dynamical systems because they generalize the notion integrating factors 
and Dulac functions. Thus, it is possible to interpret integrating factors and 
Dulac functions as invariant functions to certain cof actors: an integrating 
factor is nothing but an invariant function having cofactor 

(4.5) k = -dwX 

and a Dulac function is an invariant function with cofactor 

(4.6) k = -d[vX + j-g, 

g being a continuous function with g > for almost every z £ Q. Note 
that the latter interpretation of Dulac function as an invariant function to a 
specific cofactor makes only sense when f does not vanish, i. e., the Dulac 
function is not defined respectively singular on the invariant curve, the van- 
ishing set of the invariant function. This already gives a clue on the natural 
boundaries on the maximal domain of definition of Dulac functions. 
(2) An easy obervation is that, if f and g are invariant functions with cofactor 
k f and kg, respectively, then their pointwise product f ■ g defines also an 
invariant curve with cofactor kf + kg. 



8 



M. HIMMEL 



(3) Without loss of generality, we will always consider complex-valued invari- 
ant functions because, if f is an invariant function with cof actor k (with 
respect to some vector field in some domain), then its conjugate function f 
is also an invariant function having cof actor k, and therefore the product 
f ■ f is a real-valued invariant function with cof actor k -\-k. The same 
holds for exponential factors. 

(4) In the case of polynomial planar vector fields the algebraic part of invari- 
ant curves and exponetial factors has already been developed. We quote 
some of these results and refer to |14|, 1 17| and |[T6t for further reading. 



Therefore, let the vector field X S M^[x, y] be a polynomial vector field 
of degree d. Furthermore, assume that '^^'^^^'> -\- \ different irreducible 
invariant algebraic curves are known. Then a first integral of the form 

(4.7) H = f^'-...- 

where each fi{x, y) defines an invariant algebraic curve 

(4.8) h{x,y)=0 

for system and Aj G C not all of them null, for i = 1,2, ... ,s, 
s S N, can be constructed. The functions of type ( |4.4[ ) are called Darboux 
functions. 

(5) The irreducibilty of the invariant functions in the algebraic case must be 
replaced by the condition 

(4.9) {p G 17 : f{p) = and Vf{p) = 0} Q {p e n : X{p) = 0} 

in the non-algebraic case. 

(6) An easy obervation: Any invariant curve / = has exactly one cofactor 

^^/^- 
5. Open questions 

In section [2j proposition [1} we obtained a local existence result for Dulac func- 
tions near to a hyperbolic equilibrium. The proof was based on theorem [4j where 
we calculated a Dulac function explicitly in terms of a quadratic polynomial. In 
this proof we observed, infact, why one needs to impose that at most one eigenvalue 
of the matrix A is zero, because, if we had two such eigenvalues (7{A) = {0}, the 
matrix would have trace zero and we would have devided by it. Then, by applying 
the Hartman Grobman theorem, the result carried over to hyperbolic fixed points. 
Observe that neither does the Hartman Grobman theorem hold for equilibria p with 
linearization having purely imaginary or zero spectrum, nor could we make any use 
of it, because in both cases div X{p) = holds and our Dulac function blows up. 
Hence, in future the 

Question 2. When does a Dulac function exist near to non-hyperbolic fixed points? 

for nonlinear vector fields has to be addressed. This question is somehow naive 
because one has do deal here with the center-focus-problem, that is still not com- 
pletely solved in general. But this Dulac approach may give a new perspective on 
it. Why do we bother about all these local existence statements for Dulac func- 
tions? In principle, the motivation arose from the following 
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Algorithm 1. Input: a nonlinear planar vector field X; Ouput: number and posi- 
tion of periodic orbits together with phase portrait. 

(1) Stepl: Determine the zeros of X. 

(2) Stepl: At each zero, define locally a Dulac function. 

(3) StepS: Extend them as long as possible. If no further extension is possible, 
one has found the limit cycles and determined the phase portrait. 

Of course, this rather a pseudo algorithm because one cannot accomplish any 
of these steps. One more approachable but somehow technical step is to combine 
these local results to a global one. Technically, this means, as one has to deal 
here with different local coordinate representations, that one has to glue different 
diffeomorphisms together. 
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Fig. 1 . Gluing three different diffeomorpliisms 
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